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ESTIMATES OF. TRAPPED RADIATION ENCOUNTERED ON LOW -THRU ST 
TRAJECTORIES THROUGH THE VAN ALLEN BELTS 
by Irving M. Karp 
Lewis Research Center 

SUAAMARY 

Estimates were made of the number of trapped protons and electrons encountered by 
vehicles on low-thrust trajectories through the Van Allen belts. The estimates serve as 
a first step in assessing whether these radiations present a problem to on-board sensi- 
tive components and payload. The integrated proton spectra 4> p > E (E ranging from 
4 to 300 MeV) and electron spectra 4> > E (E ranging from 0. 5 to 4. 25 MeV) are pre- 

C 

seated for the case of a trajectory described by a vehicle with a constant -thrust accel- 
eration A equal to 10" meter per second per second. This value of acceleration cor- 
responds to a trip time of about 54 days from low Earth orbit to synchronous orbit. It is 
shown that the time spent in the belts and hence the radiation encountered vary nearly in- 

p 

versely with the value of thrust acceleration (for values of A less than about 10 

p C o 

meter/sec ). Thus, the integrated spectral values presented for the case of A = 10“° 
meter per second per second can be generalized for any other value of thrust accelera- 

2 2 o 

tion A c (less than 10 meter/sec ) by multiplying them by the factor 10 /A c . 


INTRODUCTION 

Vehicles with low -thrust electric propulsion engines are being considered for mis- 
sions from low Earth orbit to synchronous orbit or Earth escape. These spiral slowly 
through the Van Allen belts and encounter appreciable numbers of protons and electrons 
that are trapped in the belts. This radiation can be damaging to sensitive components 
and payload if it exceeds certain limits, in which case shielding has to be provided to re- 
duce this radiation. Vehicle structure and on-board equipment may provide some or all 
of the shielding required. Reference 1 discusses the damage effects of radiations from 
various sources in space and from on-board nuclear power sources on various space- 
craft subsystems. 



As a first step in assessing this problem for any vehicle, the numbers of protons 
and electrons encountered during passage through the belts have to be evaluated. Vette 
and coworkers have constructed models of the trapped radiation environment (refs. 2 
to 4) from which the proton and electron flux spectra at any position in the belts are de- 
termined. The integrated spectra of protons and electrons can then be obtained along a 
known trajectory. Reference 1 also describes calculations of such spectra for a fairly 
rapid traversal of the belts (about 1 hr) during a deep space mission. 

In this report these spectra are determined for low -thrust trajectories and general- 
ized for the case of constant thrust acceleration for all values of thrust acceleration 
o 

less than about 10 meter per second per second. The integrated proton spectrum 
4>p > E has been calculated over a range of E from 4 to 300 MeV. This extensive 
range permits accurate evaluation of the ionization dose from the protons (either un- 
shielded, or after passing through appreciable amounts of shielding). The integrated 
electron spectrum 4> e > E is calculated for energies ranging from 0. 5 to 4. 25 MeV. 

At higher energies the electron flux drops very rapidly. The electron spectral calcula- 
tions are based on a recent update of the electron environmental model (ref. 4). 


RADIATION DATA 

The model of the trapped proton environment used to evaluate the protons is de- 
scribed in references 2 and 3. These references also present orbital integration tables 
that give the average proton spectra encountered per day during circular orbits at var- 
ious altitudes (150 to 10 000 n mi) and orbital inclinations (0° to 90°). Table 3 in ref- 
erence 2, which is based on proton map AP6, gives integrated spectral data for the en- 
ergy range 4 to 30 MeV, and table 3 in reference 3, based on proton map AP7, pertains 
to the energy range from 50 to 300 MeV. 

The electrons were evaluated by using the electron environment data presented in 
reference 4. Table 9 in the reference gives the average electron spectra encountered 
per day during circular orbits at various altitudes up to 18 000 nautical miles and orbital 
inclinations from 0° to 90°. This table, based on electron maps AE4 (for the outer re- 
gion of the belts, altitudes > 10 000 n mi) and AE5 (for the inner regions), is represent- 
ative of the year 1967, a period of near maximum solar activity. During such a period 
electron fluxes in the belts are enhanced; also the fluxes in the outer regions undergo 
short-term temporal variations due to magnetic storms. Reference 4 presents time- 
averaged fluxes (averaged over several months) as well as fluxes that have probabilities 
of 0. 1 to 0. 9 of being exceeded in the puter regions. 
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TRAJECTORY 


For trajectories with low, constant tangential or circumferential thrust,* a simpli- 
fied approximate equation can be obtained (see the appendix) that presents altitude as 
a function of propulsion time for any constant value of thrust acceleration: 



where 

h altitude, meters 

R radius of Earth, 6. 380X10 6 meters 

h Q initial orbit altitude, meters 

o 

vehicle thrust acceleration, meters/sec 
t propulsion time, sec 

11 2 

p gravitational parameter , 3.986x10 cubic meters/sec 

For the purpose of evaluating radiation encountered in the Van Allen belts, this equation 

is accurate for altitudes up to synchronous orbit for values of thrust acceleration less 
-2 

than 10 meter per second per second. 

It is apparent from the equation that the time to reach a given altitude varies in- 
versely with the vehicle thrust acceleration. The total trip time and the radiation en- 
countered also vary inversely with the thrust acceleration. For the calculations, a value 

-3 

of acceleration of 10" meter per second per second (resulting in a trip time of 54 days 
from a 100-n-mi orbit to synchronous orbit) was arbitrarily selected. The radiation en- 
countered on a mission can then be generalized for any other value of acceleration A 
o 9 ^ 

(A less than 10 meter /sec ) by multiplying the calculated values by the factor 

10- S /A C . 


METHOD 

During passage through most of the trapped radiation belts the vehicle circles the 
Earth several times a day and gains a little altitude each day. Thus, the orbit integra- 
tion tables presented in references 2 to 4 give a good estimate of the radiation spectra 
encountered each day of the mission and were used to simplify the integration of 
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radiation encountered along the trajectory. For protons a set of curves was plotted 
from these tables, each curve showing the proton flux encountered per day having 
energies greater than some energy E i as a function of altitude. Curves for values of 
E i of 4, 14, 30, 50, 100, 150, 200, 250, and 300. MeV were generated. (An illustrative 
curve for cp > 14 MeV is shown in fig. 1. ) Then from the trajectory equation the alti- 
tude attained during each day of the mission was tabulated, and from each curve the flux 
encountered each day (corresponding to the daily altitudes) was determined. Table I 
shows such a listing for (p > 14 MeV and includes a running summation of the protons 
encountered during the mission. The spectra of protons encountered on both 0° and 30° 
trajectory inclinations were calculated. 

For electrons the spectra (<& e > 0. 5, 1, 2, 3, 3. 5, 4, and 4. 25 MeV) encountered on 
0° and 30° trajectory inclinations were similarly evaluated. For each trajectory two 
spectra were determined, one based on fluxes in the outer regions of the belt 
(altitudes ^ 10 000 n mi) that are time -averaged and the other based on fluxes that have a 
probability of 0. 1 of being exceeded. The latter is a more conservative estimate. Fig- 
ure 2, showing the electron flux <p e > 2 MeV encountered per day as a function of alti- 
tude, illustrates the type of data used. 


RESULTS 

Proton spectra encountered on trajectories at 0° and 30° inclinations are listed in 
table II and plotted in figure 3. Maximum numbers of protons are encountered at 0° in- 
clination. Slightly more than twice as many protons are encountered at 0° than at 30° 
inclinations. 

Above an altitude of 7000 nautical miles the fluxes of protons having energies greater 
than 4 MeV decrease sufficiently so that they contribute only slightly to the total amount 
of protons encountered on missions beyond this altitude. 

Included in table n are the ranges in aluminum corresponding to each proton energy 
listed. (The range at any energy E i is that amount of material that will stop all protons 
having energies less than E. . ) 

Electron spectra for trajectories at 0° and 30° inclinations are listed in table in and 
plotted in figure 4. Electrons encountered at 0° inclination are about 60 percent greater 
than those at 30° inclination. Also, more conservative estimates of short-term temporal 
variations in the outer regions of the belts (fluxes having probabilities equal to 0. 1 of be- 
ing exceeded) result in about 40 percent more electrons encountered than estimates from 
time -averaged fluxes. 

From figure 2 it can be seen that, although the electron flux begins to fall off with 
altitude, fairly high values of electron flux extend out to synchronous orbit altitude. Ex- 
trapolation of data beyond this orbit indicates that the flux decrease with altitude is 
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sufficient so that only a small fraction of additional electrons is encountered on missions 
beyond synchronous orbit. At synchronous orbital altitude, however, the fluxes can con- 
tribute a significant integrated flux if the period spent at this altitude is large. Table IV 
lists the electron fluxes encountered daily in synchronous orbit. These values, of 
course, are not dependent on vehicle thrust acceleration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 14, 1973, 

503-05. 
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APPENDIX - DERIVATION OF APPROXIMATE TRAJECTORY EQUATION* 


Consider the case where 

(1) The vehicle is initially in a circular orbit about the Earth. 

(2) The thrust and mass of the vehicle (and hence the thrust acceleration) are 
constant. 

(3) The thrust acceleration is small compared to the local acceleration due to grav- 
ity, so that the orbit remains approximately circular. 

(4) The thrust is directed circumferentially (it will also be nearly tangential because 
the orbit is almost circular). 

For these conditions 


torque = Fr = (1) 

dt 

where 
F thrust 

r distance from center of Earth 
m mass 
v velocity 

t time thrust has been applied 

Because the orbit is nearly circular , the centrifugal force is nearly equal to the 
gravitational force, 


mv 2 _ mq 
r ~r 2 

where p is the gravitational parameter , and 



( 2 ) 


From equation (1) 


F s A _ 1 d(vr) 
m c r dt 

' Provided by Frank J. Hrach-, NASA Lewis Research Center. 
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Substituting for v from equation (2) gives 


from which 



or 


Now 



r o = R + h o 


r = R + h 


where 

R radius of Earth 
h Q initial orbit altitude 
h orbit altitude 

Using equations (4) in equation (3) yields 



(3) 


(4) 


(5) 
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TABLE I. - CALCULATION OF INTEGRATED PROTON FLUX * > 14 MeV DURING MISSION 
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TABLE H. - PROTON SPECTRA ENCOUNTERED IN 


VAN ALLEN BELTS ON LOW-THRUST 
TRAJECTORIES 


3 2 

[Vehicle thrust acceleration A fi = 10 meter/sec ; for 
other values of A c , > Ej = (table value) 10'^/A c .] 


Orbit 

inclination, 

deg 

Energy, 

E i- 

MeV 

Proton spectra, 

*P > 'i' 2 

protons/cm 

Range (E.), a 
2 1 

g/cm of aluminum 

0 

4 

2. 9X10 12 

0. 035 


14 

1.4X10 11 

.31 


30 

2. 7xl0 10 

1.2 


50 

1.3 

2.9 


100 

7.6X10 9 

9.9 


150 

4.7 

20. 2 


200 

2.9 

33. 1 


250 

1 . 8 

48. 3 


300 

1 . 1 

65. 3 

30 

4 

1.3x10 12 

0.035 


14 

6. 0X10 10 

.31 


30 

1 . 1 

1 . 2 


50 

6. 2X10 9 

2.9 


100 

3.7 

9.9 


150 

2. 2 

20.2 


200 

1.3 

33. 1 


250 

8.2X10 8 

48.3 


300 

5.0 

65.3 


a Range (Ej) is amount of material that will completely stop 
all protons having energies less than E.. 




TABLE m. - ELECTRON SPECTRA ENCOUNTERED IN VAN ALLEN BELTS 


ON LOW -THRUST TRAJECTORIES 


o o 

[Vehicle thrust acceleration A c = 10" meter/sec ; for other values of A , 
4> e > E i = (table value) 10"®/A .] 


Orbit 

inclination, 

deg 

Energy, 

E i- 

MeV 

Electron spectra, 

O 

<f> e > Ej ; electrons/cm 

Range (E.), a 
2 L 

g/cm of aluminum 

Based on time- 
ave raged fluxes 

Based on fluxes having 
probability of 0. 1 of 
being exceeded 

0 

0.5 

2. 8x10 13 

3. 8xl0 13 

0.22 


1.0 

8. lxlO 12 

1 . 1 

.55 


2.0 

1.7 

2. lxlO 12 

1.2 


3.0 

1. 8X10 11 

2. 3X10 11 

1.9 


3. 5 

3. 1X10 10 

4. 8xl0 10 

2.2 


4.0 

3.4X10 9 

5. lx 10 9 

2.5 


4.25 

5.9X10 8 

8. 6xl0 8 

2.6 

30 

0.5 

1 . 6xl0 13 

2. 2X10 13 

0.22 


1.0 

4. 8x10 12 

6.9X10 12 

.55 


2.0 

9. 3X10 11 

1 . 2 

1.2 


3.0 

1 . 1 

1.4X10 11 

1.9 


3.5 

2.2X10 10 

3. lxlO 10 

2.2 


4.0 

2.4xl0 9 

3. 2xl0 9 

2.5 


4.25 

4.4X10 8 

5. 5X10 8 

2.6 


a Range (E.) is amount of material that will completely stop all electrons having 
energies less than E^. 



TABLE IV. - ELECTRON FLUXES ENCOUNTERED DAILY 


WHILE ORBITING AT SYNCHRONOUS ALTITUDE 


Orbit 

inclination. 


Energy, Electron flux, 

n 

Ej, cp e > Ej, electrons/(cm )(day) 


deg 

MeV 

Time -averaged 

Flux with probability of 



flux 

0. 1 of being exceeded 

0 

0.5 

2x10 11 

5X10 11 


1.0 

5X10 10 

10 11 


2.0 

3x10® 

6x10® 


3.0 

10 8 

2x10 8 


3.5 

10 7 

2X10 7 


4.0 

<10 5 

<10 5 


4.25 

~0 

~o 

30 

0.5 

9X10 10 

2X10 17 


1.0 

2X10 10 

4X 

X 

O 

h-» 

o 


2.0 

9X10 8 

2x10® 


3.0 

3xl0 7 

7X10 8 


3.5 

3xl0 6 

8xl0 6 


4.0 

© 
*— < 

V 

<10 5 


4.25 

~o 

~0 


Altitude, n mi 

Figure 1. - Proton flux <pp> 14 MeV encountered per day in circular orbit at various 
altitudes. Orbit inclination, OP. (From ref. 2, table 3. ) 

















Figure 3. - Integrated proton spectra encountered on trajectories from 
low Earth orbit to synchronous orbit. Thrust acceleration, 10" 3 
meter/ sec’ 
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Integrated electron spectra, $ e > E, electrons/cm 



(a) Orbit inclination, 0 P. (b) Orbit inclination, 30°. 


Figure 4. - Integrated electron spectra encountered on trajectories from low Earth orbit to synchronous orbit. Thrust 
acceleration, 10’’ meter/sec 2 . 
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